Identifying anomalous diffusion and melting in dusty plasmas 
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Anomalous diffusion in liquids and the solid-liquid phase transition (melting) are studied in two- 
dimensional Yukawa systems. The self-intermediate scattering function (self-ISF), calculated from 
simulation data, exhibits a temporal decay, or relaxation, with a characteristic relaxation time. This 
decay is found to be useful for distinguishing normal and anomalous diffusion in a liquid, and for 
identifying the solid-liquid phase transition. For liquids, a scaling of the relaxation time with length 
scale is found. For the solid-liquid phase transition, the shape of the self-ISF curve is found to be 
a sensitive indicator of phase. Friction has a significant effect on the timing of relaxation, but not 
the melting point. 

PACS numbers: 52.27.Lw, 52.27.Gr, 61.20.Lc, 64.70.D- 



I. INTRODUCTION 

Dusty plasma is partially ionized gas containing micro- 
sized particles of solid matter [l[. In a plasma the 
sheath above a lower electrode has electric fields that 
can levitate and confine highly charged particles, so that 
they are suspended. When only a single layer is sus- 
pended, the interaction between dust particles is a repul- 
sive Yukawa potential ■ Video microscopy allows imag- 
ing this two-dimensional (2D) suspension at an atomistic 
scale, so that we can track particles and measure their 
individual positions and velocities in each video frame, 
yielding the same kind of data as the molecular dynam- 
ics (MD) simulations, reported here. 

Particles self-organize in a crystal-like triangular lat- 
tice with hexagonal symmetry due to strong interpar- 
ticle interaction. In this strongly-coupled plasma, the 
Coulomb interaction with nearest neighbors is so strong 
that particles do not easily move past one another [3j . In 
many experiments, the particles occupy only a single hor- 
izontal layer, and are often described as 2D experiments 
P, US] • Dusty plasma is a driven-dissipative system [lj , 
and its kinetic energy is determined by the balance of the 
energy input and dissipation. As the driven energy input 
increases, the lattice becomes disordered in a solid-liquid 
phase transition 0, [H-Qjl . 

Random particle motion in dusty plasmas can be di- 
vided into several stages. Ballistic motion 14[ occurs on 
a short time scale < while caging oscillations [HI 

happen on a typical time scale of 10 L)~ d . Here uj v d is the 
nominal plasma frequency [lit ]. At later times, particles 
can escape their cages and diffuse. A current research 
topic that has not been resolved is whether this long- 
time random motion exhibits normal diffusion or anoma- 
lous diffusion. Experiments [3, 12 1 and simulations 
lH-2o[ with these 2D systems have indicated superdif- 
fusion, where the mean-square-displacement (MSD) in- 
creases with the time more rapidly than linear scaling, 
but other simulations suggest that random motion may 
be normal diffusion at sufficiently long times [2l|. Iden- 



tifications of normal diffusion and anomalous diffusion 
(e.g. superdiffusion here) have typically made use of dy- 
namical measures: MSD time series, the probability dis- 
tribution function (PDF) for particle displacements, and 
the velocity autocorrelation function (VACF) [18l42l| . 

Data analysis methods used to study the solid-liquid 
phase transitions can be grouped in two categories, static 
and dynamic. Here, we term a method static if the input 
data can be a single snapshot of particle positions, but 
dynamic if it requires a series of positions for each par- 
ticle. According to this classification, most attempts to 
identify phase transitions in dusty plasma experiments 
and MD simulations have employed only static struc- 
tural measures, such as Voronoi diagrams and correlation 
functions for particle positions and angular orientation 
0, IH-EH [22| ■ Sheridan used a static method of applying 
the empirical Lindemann criterion 13 , 23[ . A dynamical 
method of identifying melting has been developed theo- 
retically, using empirical criteria based on the long-time 



and short-time self-diffusion coefficients 24. 251; this has 



been demonstrated to be useful for simulations of both 
2D [H and 3D [H] systems. 

Here we carry out dynamical analysis using the inter- 
mediate scattering function to characterize random mo- 
tion. We study two physical processes: anomalous diffu- 
sion in liquids and the solid-liquid phase transition. 

In Sec. II, we will review the intermediate scattering 
function briefly. In Sec. Ill, we will introduce our two MD 
Yukawa simulation methods: Langevin and frictionless. 
They will model a 2D dusty plasma. In Sec. IV, we will 
present results for the two physical processes: normal 
and anomalous diffusion in liquids, and the melting phase 
transition. 



II. SELF-INTERMEDIATE SCATTERING 
FUNCTION 

The intermediate scattering function (ISF) (2(| , which 
has been used widely in other fields, is also called the 
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density-density correlation function [27|. The ISF is de- 
fined in terms of the particle trajectories: 

f ^ 4^^ (exp( " ik ' [ri(t) ' ri(o)l)) - (i) 

i i 

Here, Vi(t) is the trajectory of the ith particle in the 
system consisting of N particles. The Fourier transform 
variable k is usually called a wave number, although no 
waves are studied using this method. Equation (JTJ) makes 
use of an ensemble average (•■■), which in practice is done 
by averaging for various initial starting times in place of 
t = 0. 

Calculating the ISF directly from Eq. (JTJ) is a dynami- 
cal analysis method because it requires as its input data 
a time series of positions for each particle. This is the 
method we will use, as has been done previously in ex- 
periments with granular materials [28|, [2^| and in MD 
simulations [30l- 32 1 of systems other than dusty plasmas. 

Besides starting from measurements of particle trajec- 
tories, other experimental methods of obtaining the ISF 
have been devised for colloids [HI, supercooled liquids 
34], and polymer nanocomposites [HI . In these exper- 
iments, the ISF was determined from data produced by 
dynamic light scattering [33j, x-ray photon correlation 
spectroscopy 35, 3(|, or neutron spin echo spectroscopy 
[34I I37I ] . A comparable experiment with a dusty plasma 
was reported by Khodataev et al. [38| , using photon cor- 
relation spectroscopy to yield a function of wave number 
and time. 

It is useful to compare the ISF, which is a time series, 
with two related functions. First, the dynamic structure 
factor [37] is the Fourier transform of the ISF; it is defined 
as Sik, to) = f_ 0O F(k,t) exp(iwi) dt/2it. This dynamic 
structure factor has been used for studying 2D Yukawa 
liquids in both theory [3j| and simulations [40j. Second, 
phonon spectra, which are widely used for experimental 
and simulation data, have two parts, for longitudinal and 
transverse waves. The longitudinal phonon spectrum is 
related to the dynamic structure factor [4l[, although 
phonon spectra are more commonly computed differently, 
as a Fourier transform of the particle current [42j . 

The ISF is composed of two parts [13, S3], F(k,t) = 
F s (k,t) + F c (k,t). The most commonly used part is 
F s (k,t), which is often called the incoherent part, or the 
self-ISF: 

F s (k,t) = l^( ex p(-ik. [ ri (t) -ri(0)])). (2) 

i 

The less commonly used part is F c (k,t), which is called 
the coherent part: 

F c (k,t) = ^£X>x P (-ik- Kt) - rj (o)])>. (3) 

The self-ISF, F s (k, t), is a measure of single-particle 
dynamics as a function of time. This makes it compara- 
ble to the MSD and PDF, which are also computed from 



the trajectories of individual particles recorded for a long 
time. Thus, the self-ISF can be used to study some of 
the same physical phenomena as MSD and PDF, such as 
random motion and the related idea of relaxation [29| . If 
random motion consists of normal diffusion, as for exam- 
ple with Brownian motion with a diffusion coefficient D, 
then [H, [H, S3| 



F s (k,t) ~ exp(-Dfc 2 i). 



(4) 



In Sec. IV A 2, we will generalize Eq. (H|) for the case of 
anomalous diffusion, such as superdiffusion. 

The self-ISF, Eq. J2]), is often used by itself, without 
reporting the coherent part, Eq. ([3]). This is a common 
practice with particle trajectory data from experiments 
|28l [29| and MD simulations [3CH32| for various physical 
systems. Here, we will also use only the self-ISF for our 
Yukawa simulations of dusty plasmas. 

A graph of the self-ISF typically reveals two stages of 
random motion. We illustrate this in Fig. 1 with a sketch 
of F s (k, t) for a normal liquid (not supercooled). The 
decay of the F s (k, t) curve is the signature of what is often 
called relaxation. Caging motion is indicated at short 
times. This early part of the curve is sometimes termed 
the fast (3 relaxation. Diffusion is indicated at long times, 
as particles gradually escape their cages [29| . In this later 
part of the curve, sometimes termed a relaxation, F s (k, t) 
gradually decays toward zero. Here, we will use the self- 
ISF two ways. We will inspect the shape of the self-ISF 
decay, and sometimes fit it to a stretched exponential 
44]:" 



F s (k,t) = exp[-(Vr(fc))^ fc )], 



(5) 



where r(fc) is a relaxation time. (Note that this use of 
the symbol (3 has no relation to the /3 relaxation.) 



III. SIMULATION 



A. Parameters 



Equilibrium Yukawa systems are characterized by two 
dimensionless parameters: the coupling parameter F 
and the screening parameter k 45j, |46j. Here, V — 
Q 2 /(47reoafcsT) and k = a/Xjj, where Q is the parti- 
cle charge, T is the particle kinetic temperature, Ad is 
the screening length, a = (nir)^ 1 is the Wigner-Seitz ra- 
dius 16j, and n is the areal number density. Another 



characteristic length is the lattice constant b for a defect- 
free crystal, which is b — 1.9046 a for a 2D triangular 
lattice. 

Our two simulation methods are the same in many 
respects. Both simulation methods use a binary inter- 
particle interaction with a Yukawa pair potential, 



<f>i,j = Q 2 (47re r lJ ) 1 exp(-r i j / X D ), 



(G) 
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where r^j is the distance between the ith and jth par- tions [l8|. The equation of motion is 
ticles. In both simulations, particles are only allowed to 

move in a single 2D plane. Conditions remained steady mr^ = ~V / J^jji 

during each simulation run. For both simulations, the 
parameters we used were TV — 16 384 particles in a rect- 
angular box with periodic boundary conditions. The box 
had sides 137.5 6 x 119.16. We truncated the Yukawa 
potential at radii beyond 12 6 [47]. The integration time 
step was 0.037 and simulation data were recorded 
for a time duration of 1777 lu~^ in after the system 
reached its steady state. Other simulation details are 
presented in H, 47|. We report results with distances 
normalized by 6, while time (and frictional damping rate 
v) are normalized using the nominal plasma frequency 
u p d = (Q 2 /27re wa 3 ) 1 ^ 2 [3], where m is the particle 
mass. 

We will next review the two simulation methods. They 
differ mainly in the equations of motion that are solved. 



(8) 



B. Langevin MD simulations 



which we integrate for all particles. A Nose-Hoover ther- 
mostat is applied to maintain a desired temperature [l8| . 

This MD simulation method describes a frictionless 
atomic system. The particles collide among themselves, 
without any interaction with gas or other external influ- 
ences. This method mimics thermal equilibrium condi- 
tions. 

There are two parameters we can change in the fric- 
tionless MD simulations: L and k. In the Langevin sim- 
ulations, we can also vary v. Varying T and k is equiva- 
lent to varying temperature and density, and we will vary 
them over a range that allows us to simulate liquids or 
solids. 

Our method is to generate trajectories i"j(t) for all par- 
ticles by integrating Eq. ([7]) or ([5]) , and then to compute 
the self-ISF using Eq. ([2]). The self-ISF is a time series. 
We repeat its calculation for various wave numbers, k. 



Our Langevin MD simulations take into account the 
dissipation due to frictional gas damping. The Langevin 



equation [18M21L 1471 . |48j of motion for each particle is 



mr, 



= -v£- 



(7) 



Trajectories rj(t) are found by integrating Eq. ([7]) for 
all particles. Terms on the right-hand side include a 
frictional drag vmii and a random force d(t). Note 
that we retain the inertial term on the left-hand-side in 
Eq. (0, unlike some Brownian-dynamics simulations of 
overdamped colloidal suspensions [lij], where it is set to 
zero. 

Our Langevin simulations mimic 2D dusty plasma 
experiments [1], but the driven-dissipation mechanism 
is only an approximation of the processes in experi- 
In our Langevin simulation, the heating 



47]. 



ments 

and friction are explicitly coupled by the fluctuation- 
dissipation theorem l51| : this models collisions with 
gas atoms that provide both frictional drag and random 
kicks. However, besides random kicks from gas atoms, in 
dusty plasma there are some additional heating mech- 
anisms arising from ion flow, particle charge fluctua- 
tio ns [52l . l53| , and sometimes external laser manipulation 
[3, l47| that are not explicitly modeled in our Langevin 
simulations. 



IV. RESULTS AND DISCUSSIONS 

We present results for two physical processes for 2D 
systems under steady conditions. First, we evaluate 
whether random motion in a 2D liquid exhibits normal 
or anomalous diffusion. Second, we test whether the self- 
ISF can serve as a sensitive indicator of the phase tran- 
sition between solid and liquid. 

In this paper, relaxation refers to the decay of the self- 
ISF, F s (k, t), as shown in Fig. 1. The relaxation rate 
in general depends on the scale length, which is param- 
eterized here by 1/k. U sing the terminology of other 



users of the self-ISF [29|, |54|, the early stage of decay 



is termed ft relaxation; during this early time particles 
are mainly trapped within their cages formed by nearest 
neighbors. The later stage is termed a relaxation, and 
this corresponds to diffusion as particles decage. The 
term "decaging" refers to a particle's movement so that 
it is no longer trapped by the previous nearest neighbors. 
In a liquid, particles decage much more rapidly than in 
a solid, so that this ft relaxation is much faster in liquids 
than solids. 



A. Anomalous diffusion in 2D liquids 

1. Results for the self-ISF 



C. Frictionless equilibrium MD simulations 

In addition to our Langevin MD simulations which in- 
clude friction, to obtain results in the frictionless limit, 
we also performed frictionless equilibrium MD simula- 



Results from our Langevin simulations are presented 
in Fig. 2 for typical liquid conditions far from the phase 
transition. Curves are shown as functions of time for 
various values of k. Note the smooth and gradual decay 
from unity to zero as time increases. This decay is the 
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signature of relaxation. Here, the decay develops without 
any plateau. The lack of a plateau in the time-variation 
of the self-ISF is similar to what is seen in granular flows 
[28l [2^], but different from what is expected for super- 
cooled liquids and glasses (cf. Fig. 3 of (H3|). 

To help quantify the relaxation that is observed in 
Fig. 2, we fit the time-dependence of the self-ISF to 
the empirical form Eq. ([5]). Since the relaxation process 
spans many decades of time, to perform this fit without 
biasing results toward long times, data points were sam- 
pled from the simulation results at time intervals equally 
spaced on a logarithmic scale. For a liquid far from the 
phase transition, Eq. ([5]) fits our simulation data points 
well (shown as solid lines in Fig. 2). The two free pa- 
rameters for the fit, r(fc) and /3(k), help quantify the 
relaxation process. We discuss their physical significance 
below. 



2. Searching for anomalous diffusion 

We develop two tools for identifying anomalous diffu- 
sion. Previous investigators have usually used the MSD, 
looking for a scaling with time that differs from the MSD 
oc t scaling expected for normal diffusion 3, 12 -III 21 1. 
For 2D systems, data are typically fit to the form 



(r 2 (t)) = 4D t a , 



(9) 



where a = 1 is the case of normal diffusion, a > 1 is 
superdiffusion, and a < 1 is subdiffusion. (This use of the 
symbol a has no relation to the a relaxation mentioned 
above.) Here we introduce two other tools that are also 
based on how random motion develops with time: the 
scaling with r vs. k, and the value of (3. Recall that r 
and j3 are the fitting parameters for Eq. (0. 

Our first new tool is the power-law scaling of the fitting 
parameter r as compared to fc. To do this, we must first 
generalize Eq. §4$ to allow for anomalous diffusion. Start- 
ing from Eq. ([2]) , previous authors [43l . |55| have demon- 
strated that 



F s (k,t)~exp( V^), 



(10) 



where we have substituted 4 in place of 6 in the denomi- 
nator for two dimensions instead of three. Next, we sub- 
stitute Eq. © in Eq. (flQ]t . yielding 

F s Ql, t) ~ exp(~k 2 Dt a ) = exp(-D(k 2/a t) a ). (11) 

Examining the argument on the right-hand-side reveals 
the scaling 



T OC k- 2 ' C \ 

In the case of normal diffusion, a 
t oc fc -2 , as previous authors have noted 28|, |29j, |4 



(12) 

1, the scaling is 



Here, we note that the superdiffusion case a > 1 has 
r varying with a lesser power. Thus, the signature of 
superdiffusion will be a slope weaker than —2 when r 
is plotted vs. k using log-log axes. Similarly, a slope 
stronger than —2 in the same r — k plot is the signature 
of subdiffusion. 

Our second new tool is the fitting parameter (3 for the 
sclf-ISF. Comparing Eq. ([5]) and Eq. (iTl) , we see that 
the value of (3 is essentially the same as a. The only 
difference is that when using actual data, the value of (3 
is generated by a fit, while a is generated by examining 
a log- log plot. The community of scientists who use the 
sclf-ISF traditionally uses (3, although until now it has 
not been used as an indicator of superdiffusion. Scientists 
who use MSD to characterize superdiffusion, on the other 
hand, traditionally use a. 

The two new diagnostic tools we described above may 
be useful for experiments in other fields where measure- 
ments allow one to obtain the self-ISF. While methods 
exist to distinguish anomalous diffusion from normal dif- 
fusion (using MSD, PDF, or VACF) if particle track- 
ing is feasible, as in dusty plasma experiments 1J| and 
MD simulations 18l421j , other methods are needed when 
tracking is impossible. For example in a simple liquid, the 
motion of individual molecules cannot be tracked, but 
the self-ISF can be obtained using dynamic light scat- 
tering [43| or some other spectroscopic methods [34T4371 ] . 
Thus, the two diagnostic tools presented above could find 
an application in such experiments. 



3. Results for fitting the self-ISF 

For the conditions of a liquid far from the phase tran- 
sition, we use our two tools (scaling of r vs. k and the 
value of /3) to test for anomalous diffusion. We use our 
two MD simulations, Langevin and frictionless. Ott and 
Bonitz [21| previously varied the values of friction v and 
observation time over wide ranges, and using the MSD 
method prepared a diagram showing the conditions that 
favor normal diffusion or superdiffusion. This diagram, 
Fig. 3 of [2l|, predicts that the value of v which we use in 
our Langevin simulation will yield normal diffusion, while 
the v = case of our frictionless simulation will yield su- 
perdiffusion over any reasonable observation time. Here 
we test whether our two cases, analyzed using our two 
new tools, yield the same conclusion as in [21] . 

For the r vs. k scaling, in Fig. 3(a) we find the scal- 
ing r oc k~ 2 for the Langevin MD simulations, and the 
scaling of r oc fc -7 (7 < 2) for the frictionless MD sim- 
ulations. In other words, random motion is diffusive for 
our frictional (Langevin) case, but superdiffusive for our 
frictionless case. This is in quantitative agreement with 
Fig. 3 of [2l|], prepared using MSD curves, which demon- 
strated that friciton can inhibit superdiffusion 1|| 21 1. 

For /3, in Fig. 3(b) we find values near unity for our 
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Langevin simulation, but a value definitely > 1 for the 
frictionless simulation, for moderate values of k. This is 
again consistent with the conclusion of diffusive motion 
for our frictional (Langevin) case, but supcrdiffusivc mo- 
tion for the frictionless case. At extremely small or large 
values of k, however, (3 can be different. Figure 3(b) re- 
veals an overall trend for /3 to increase with k, especially 
at extremely small or large values of k. Previous authors 
28, 29] have identified dynamic heterogeneities as the 
cause for (3 < 1 for very small k, i.e., very large length 
scales. For large k, previous authors have not reported 
enhanced values of (i like those we see in Fig. 3(b). One 
possible interpretation of our large k observation is that, 
at these short length scales, the self-ISF is affected more 
by caging motion than by random walks associated with 
decaging. 

To summarize, we find that the relaxation of the self- 
ISF is a sensitive indicator to distinguish anomalous dif- 
fusion from normal diffusion. The indication can be made 
using either of the two fitting parameters, t and j3. 



B. Solid-liquid phase transition 

Simulation studies of the solid-liquid phase transition, 
which is also called an order-disorder transition or melt- 
ing, are generally done using measures of static structural 
order, such as defects or correlation functions of particle 
position or bond orientation [H, H|-13, 22 1. Dynamical 
measures can provide additional information that can be 
helpful for identifying a phase. Temperature is kind of 
dynamical measure for random motion, and so is the self- 
ISF. 

Our first goal, for phase transitions, is to test the use 
of the self-ISF as an indicator of the phase transition. 
We perform tests that indicate that it is sensitive in dis- 
tinguishing solids and liquids near the phase transition. 
This development is useful because the self-ISF is based 
on dynamics rather than structure. Our second goal is 
to determine what role friction plays in the phase tran- 
sition. We will vary temperature and density, using the 
normalized quantities T and k, and we will also vary the 
friction v, to determine whether the self-ISF is sensitive 
to phases, and what role friction plays. 



1. Dependence onY 



We note two features of the self-ISF curves in Fig. 4 
that are different on either side of this sudden change. 
First, the gap between curves is much wider for low tem- 
perature (high r) conditions in the upper right of the 
figure as compared to the high temperature conditions in 
the lower left. We varied T in small steps near Y = 200, 
where we find the sudden change in the gaps between the 
self-ISF curves. Second, we identify different shapes for 
the decay of the self-ISF for low and high temperatures. 
For the high temperature (low Y) conditions expected for 
liquids we found, in Sec. IV A, that the self-ISF decays 
according to the empirical law Eq. ([5]), but for low tem- 
peratures we found that Eq. ((5]) does not come even close 
to the shape of the curves in the upper right of Fig. 4. 

Comparing to previous simulations that used struc- 
tural measures, we can confirm that the sudden change 
in the self-ISF curve corresponds to the phase transition. 
Using a measure of local orientation order that exhibited 
a large jump at the phase transition, a phase transition 
curve for Y vs. k was reported, Fig. 6 in 22[. Inter- 
polating their results, we find that the phase transition 
occurs at about Y = 200 for k = 1.2, for a 2D Yukawa 
system modeled with a frictionless MD simulation. This 
result is consistent with the sudden change that we ob- 
served in our self-ISF curves for k = 1.2 in Fig. 4: first, 
the curve's shape is changed; second, the curves are nar- 
rowly spaced for liquids (r < 200) and widely spaced for 
solids (r > 200). (One difference in the simulations of 
[22| and ours is our use of friction comparable to values 
in 2D dusty plasma experiments. We will explore the role 
of friction in Sec. IV B 4.) 

Thus, we conclude that the self-ISF curve is very sen- 
sitive to phase. It shows promise to become a reliable 
indicator of the phase transition, although further tests, 
for different physical parameters would be needed to con- 
firm its reliability. 

An attraction of this method is that the self-ISF is a 
dynamical rather than structural measure. Combining 
both dynamical and structural measures can be useful in 
distinguishing phases. For example, supercooled liquids 
have the structure of normal liquids, so that a dynam- 
ical measure is needed to distinguish the two. We note 
that other methods of using dynamical measures for such 
purposes include the Lowen criterion 24| 25|, which in- 
volves a comparison of diffusion coefficients computed in 
two limiting cases. 



We find that the self-ISF time series undergoes a sud- 
den change at the melting point. This is seen for the 
self-ISF in Fig. 4 for our Langevin simulation. Here, we 
have chosen to present results for a small wave number 
k = 27r/6, which is the wave number corresponding to a 
lattice constant so that the self-ISF indicates dynamics at 
the length scale of nearest neighbors. As we varied Y in 
Fig. 4, we held k — 1.2 and v/uj p d = 0.027 as constants. 



2. Dependence on n 

Varying only the density or k, we again find the same 
two results as for varying Y. There is a sudden change in 
the gaps between curves, and the curves take a different 
shape at a point that we can identify as the phase tran- 
sition. This is seen in Fig. 5, where we changed k in our 
Langevin simulations with Y = 200 and v/uj p d — 0.027. 
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The transition occurs at n 
Fig. 6 in p2| . 



1.2, which is consistent with 



3. Structure relaxation time 

We investigate how rapidly disorder develops on the 
length scale of a cage, i.e., the interparticle spacing b. 
We will test whether it occurs with an Arrhenius or 
Vogel-Fulcher law as in other complex fluids such as col- 
loids [5? !, g ranular materials [29| , supercooled liquids and 
glasses [57|] . We carry out this investigation by charac- 
terizing a decay time for the self-ISF. We could use the 
fit parameter r(k) from Eq. (O, but for simplicity, rather 
than fitting the self-ISF, here we will adopt the practice 
of other authors of measuring the time required for the 
self-ISF to decay by a factor of 1/e [56. 57|. This is done 
only for k — 2n/b, corresponding to the length scale of 
a cage. Following the terminology used for example in 
the literature for colloids [56j |. supercooled liquids and 
glasses [13], this 1/e decay time, denoted here as T 2ir /b, is 
called the "structure relaxation time." It is in principle 
the same as r if j3 = 1. 

Results for the structure relaxation time r 2Tr /b for the 
curves of Figs. 4-5 are shown in Fig. 6, revealing how 
T 27r/f) varies with temperature and number density. As 
the normalized temperature 1/r increases, the system 
melts and r 2ff /6 decreases about one order of magnitude. 
The rate at which a structure relaxes, l/T 2ir /b, increases 
with temperature, which is the same trend, for exam- 
ple, as for diffusion 58]. Melting also occurs as k is in- 
creased. Plotting T 27r /b vs. number density, in the inset 
of Fig. 6(b), we note a nearly linear scaling in the liquid 
regime. 

The scaling we observe for T 27r /b vs. number density 
is different from the Arrhenius and Vogel-Fulcher laws. 
The Vogel-Fulcher law, which has been found empirically 
for other complex fluids 29|, 56, 57], has r 27r / h diverging 
to infinity as the number density increases toward the 
phase transition point. Here, T 27r / fc increases about an 
order of magnitude near the phase transition, but it does 
not diverge. 



4- Dependence on v 

We find that friction can play an important role in 
the random motion of 2D Yukawa liquids, but only at 
large damping rates. This result is shown in Fig. 7 where 
we varied v while holding T and k constant in a liquid 
regime. We observe that the self-ISF curves all lie on 
top of one another for v/io p d J$ 0.03, but not at higher 
damping rates where they decay more slowly as v is in- 
creased. The curves at higher damping rates have the 
same general shape, but are retarded in time. In other 
words, friction has a significant effect on the timing of re- 



laxation, but not the melting point. For the limiting case 
of no friction, we also include in Fig. 7 results for our fric- 
tionless equilibrium MD simulation, and these agree with 
the Langevin simulation for low friction v/uj p( i < 0.03. 
We note that Vaulina et al. [58j found that diffusion is 
mostly independent of damping rate, in the limit of low 
friction, which is similar to our results for relaxation. 

We interpret the difference in the F s (2ir/b, t) curves at 
various damping rates as indicating retardation of diffu- 
sion at high friction levels. At higher damping rate, the 
random motion of particles is resisted, then more energy 
is dissipated locally. As a result, the collective relaxation, 
which refers to the diffusion, will be retarded. 



V. CONCLUSIONS 

In conclusion, we have performed Yukawa MD simula- 
tions to study two physical processes in 2D dusty plas- 
mas: anomalous diffusion in liquids and melting. For 
liquids, examining the decay or relaxation of the self- 
ISF reveals a scaling of the relaxation time vs. length 
scale. This scaling is demonstrated to be useful for dis- 
tinguishing normal and anomalous diffusion. The self- 
ISF is also demonstrated to be a sensitive indicator of 
the solid-liquid phase transition, i.e., melting. Friction 
has the effect of retarding relaxation. 

This work was supported by NSF and NASA. 
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FIG. 1: A typical example of the self-intermediate scatter- 
ing function (self-ISF) F s (k,t), sketched here to show how 
it begins at unity at t = 0, and then relaxes gradually to 
zero as t — Y oo. Relaxation happens in two steps: early-time 
caging motion, and long-time diffusion. The self-ISF curve is 
often modeled by Eq. (0). The case shown here is for a nor- 
mal liquid in 2D, and to reveal phenomena at a length scale 
corresponding to caging we chose k = 3.3 a, where a is the 
Wigner-Seitz radius, as defined in Sec. III. 

FIG. 2: (color online). Time dependence of the self-ISF for 
various wave numbers (fc) for the Langevin molecular dynam- 
ics (MD) simulation in the liquid regime: V — 200, n — 2 
and v/ujpd = 0.027. The solid lines are the corresponding fits 
to Eq. ©. Quantities are normalized by parameters, such as 
the lattice constant 6, as defined in Sec. Ill A. 

FIG. 3: (color online). The fitting parameters, relaxation 
time r in (a) and exponent j3 in (b), as a function of k. 
These cross symbols are from fitting the data in Fig. 2 for 
the Langevin MD simulation, and the circle symbols are for 
the frictionless MD simulation for the same liquid conditions, 
r = 200, k — 2. We find that the scaling in (a) and the values 
of /3 in (b) are useful as indicators of anomalous diffusion. 

FIG. 4: (color online). The self-ISF for the length scale cor- 
responding to the lattice constant, k = 2n/b. Results shown 
are for the Langevin MD simulations with constant values of 
K — 1.2 and v/u} p d = 0.027. To search for an indication of a 
phase transition, T is varied over a range spanning both the 
solid and liquid phases. A phase transition near Y — 200 is 
indicated by a change in the curve's shape, and by larger gaps 
between curves for solids, T > 200. 

FIG. 5: (color online). The self-ISF at k = 2ir/b, as in Fig. 4 
but with varying k. The simulations here are all Langevin 
MD simulations with constant values of F = 200 and v /tj p d = 
0.027. To search for an indication of a phase transition, k is 
varied over a range spanning both the solid and liquid phases. 
A phase transition near « = 1.2 is indicated by a change in 
the curve's shape, and the gaps between curves. 

FIG. 6: Structure relaxation time measured from Figs. [4] and 
[S]as a function of (a) temperature and (b) screening parameter 
k. The structure relaxation time in our 2D underdamped 
Yukawa system does not obey the Arrhenius law or Vogel- 
Fulcher law. The same data are plotted in the inset as a 
function of number density 1/k 2 , for the liquid regime. 

FIG. 7: (color online). The role of friction in relaxation is in- 
vestigated by computing the self-ISF at k = 2n/b, for simula- 
tions of a liquid. The simulations here include both Langevin 
MD simulations and frictionless MD simulations in a liquid 
regime, with constant values of T = 200 and n — 3. Unlike 
temperature and number density, changing the friction does 
not significant change the shape of the curves. The effect of 
friction in the self-ISF curves is a retardation of the decay at 
higher damping rate. 
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